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Dear  Dr.  Rood: 


Free-Surface  Effects  on  Shin  Boundary  Lavers  and  Wakes 
ONR  Contract  NOOQU-SS-K.QUS 
Final  Technical  Report 

This  is  the  Final  Technical  Report  on  the  above  referenced  research  project  through 
31  December  1991.  Future  related  research  will  be  reported  under  the  recently  initiated 
contract  renewal  under  ONR  Grant  N00014-92-J-1092. 

1.  Foil-Plate  Model  Experiments 

The  report  and  journal  article  describing  the  wake  experiments,  including 
comparisons  with  the  calculations,  has  been  completed  (Stem  et  al.,  1991). 

— A  new  experimental  program  has  been  planned,  which  will  focus  on  the  detailed 
resolution  and  determination  of  the  fundamental  physical  mechanisms  of  the  flow  in  the 
neighborhood  of  the  body/free-surface  boundary  layer  and  wake  (i.e.,  the  solid-fluid 
juncture  and  wake  regions).  Both  mean  and  turbulence  velocities  will  be  measured  for 
zero,  medium,  and  large  steepness.  Also,  an  effort  will  be  made  to  measure  the  wave- 
induced  separation  flow  pattem>,^e  experiments  require  a  substantial  upgrade  of  the  IIHR 
towing  tank,  including  the  folloWtg:  two-component  LDV  system  designed  specifically 
for  the  present  project,  which  utilizes  t  streamlined  fiber-optic  pri  be  and  an  automated 
traversing  system;  complete  rebuilding  of  the  towing-tank  trailer  in  order  to  house  the  LDV 
system;  rewiring  of  the  towing  tank  to  provide  the  necessary  power  requirements  of  the 
li)V  system;  and  complete  refurbishing  of  the  foil-plate  model.  The  tot^  cost  involved  is 
about  $200  K.  Presendy,  most  of  the  issues  pertaining  to  the  upgrade  have  been  resolved. 
We  hope  to  complete  the  upgrade  and  initiate  the  experiments  in  about  six  months. 

2.  Series  60  C£  =  .6  Shin  Model  Experiments 

The  report  and  journal  article  describing  the  mean-flow  measurements  in  the 
boundary  layer  and  wake  and  wave  field  for  Ft  =  .16  and  .316  has  been  completed  (Toda 
et  al.,  1991).  The  journal  article  describing  the  scale-effect  studies  has  also  been  completed 
(Longo  et  al.,  19S>91). 
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The  local  wave-elevation  measurements  have  been  completed  for  5  and  10  deg  yaw 
angle  and  Fr  =  .16  and  .316.  Typical  results  are  shown  in  figure  1,  which  are  considered 
satisfactory  and  dramatically  underscore  the  complex  nature  of  the  wave  field  for  the  yawed 
condition.  Presently,  comparisons  are  being  made  between  the  data  and  inviscid-flow 
calculations.  Also,  the  det^ed  mean-velocity  and  pressure  field  measurements  are  being 
initiated  for  the  10  deg  yaw  condition. 

3.  Development  of  the  Compufationai  Method 

During  the  present  contract  period,  our  efforts  have  been  directed  in  the  following 
areas:  (a)  calculations  for  the  Stokes- wave/flat-plate  flow  field;  (b)  development  of 
interactive  and  large-domain  approaches  for  calculating  ship  boundary  layers  and  wakes  for 
nonzero  Fr,  and  (c)  extensions  of  the  interactive  approach  for  yawed-ship  calculations.  The 
results  from  these  efibrts  are  summarized  below. 

(a)  Laminar-flow  calculations  are  in  progress  concerning  higher-order  and  exact 
treatments  of  the  free-surface  boundary  conditions  (fsbc)  and  detailed  resolution  of  the  flow 
in  the  region  very  close  to  the  free  surface  in  which  the  fsbc  play  an  important  role.  Results 
for  wave  steepness  Ak  =  .01  have  been  obtained  utilizing  a  refined  grid  near  the  free 
surface  and  three  alternative  treatments  of  the  fsbc:  (1)  zero-gradient  conditions  for  both  V 
and  p;  (2)  zero-gradient  conditions  for  V  and  an  inviscid  small-amplitude  wave  normal- 
stress  condition  for  p;  and  (3)  exact,  i.e.,  viscous  stress  conditions  for  u,  v,  and  p, 
continuity  for  w,  and  free-surface  update  through  the  kinematic  fsbc. 

Figure  2  provides  a  comparison  of  the  relative  error  for  the  three  treatments, 
including  results  in  which  the  conditions  are  satisfied  both  on  the  mean  and  updated  free 
surfaces.  Clearly,  the  error  is  substantially  reduced  for  the  more  exact  treatments.  Figures 
3a, b  provide  an  overview  of  certain  aspects  of  the  solutions  through  the  velocity  [AV  = 
V(Ak  *  .01)  -  V(Ak  =  0)]  and  vorticity  [A<a  =  a)(Ak  =  .01)  -  (ofAk  =  0)]  differences  for 
zero  and  nonzero  steepness.  The  effects  of  the  wave-induced  pressure  gradients  are  seen  to 
propagate  to  depths  of  about  Xll  and  primarily  influence  Au,  Aw,  and  Acox.  and  Acoz-  The 
effects  of  the  fsbc  are  confined  to  a  region  very  close  to  the  free  surface  and  primarily 
influence  Au  ,  Aw,  and  Ao).  In  particular,  the  fsbc  induce  vorticity  and  vorticity 
interactions,  as  displayed  in  figures  4a,b,c,  which  display  the  convection, 
stretching/bending,  and  diffusion  terms  of  the  axial  tOx,  transverse  cOy,  and  depthwise  coz 
vorticity  equations,  respectively. 

Presently,  the  solutions  are  being  further  analyzed  to  fully  explicate  the  role  of  the 
fsbc  and  the  interaction  between  the  boundary-layer  and  wake  vorticity,  including  wave- 
induced  pressure  gradient  effects,  and  the  vorticity  induced  by  the  fsbc  and  the  free-surface 
boundary  layer  edge  conditions.  Also,  calculations  using  exaggerated  free-surface 
boundary-layer  edge  conditions  and  for  larger  wave  steepness  (i.e.,  Ak  =  .11)  are  in 
progress.  An  abstract  for  a  proposed  presentation  describing  this  work  has  been  submitted 
to  the  1992  Free-Surface  Vorticity  Workshop  (see  item  4). 

(b)  Validation  studies  have  been  completed  for  the  interactive  approach  through 
calculations  for  the  Series  60  Cb  =  .6  ship  model,  including  resolution  of  the  bow  flow 
and  comparisons  with  the  experimental  data  of  Toda  et  al.  (1991).  The  results  are 
satisfactory  and  indicate  that  the  interactive  approach  fully  exploits  the  advantages  of  both 
the  inviscid  and  viscous  technologies  and  properly  accounts  for  their  interaction. 
Presently,  a  journal  article  is  in  preparation  describing  the  calculations  and  comparisons 
(Tahara  and  Stem,  1992). 
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With  regard  to  the  development  of  a  large-domain  approach  for  the  present  problem 
two  strategies  for  initiating  the  solution  have  been  pursued:  (1)  based  on  the  nonzero  Fr 
inviscid  solution;  and  (2)  based  on  the  zero  Fr  viscous  solution.  As  discussed  in  the  last 
Quarterly  Progress  Report  dated  19  July  1991,  satisfactory  results  were  obtained  using 
strategy  Q)  and  both  fixed  and  free-surface  conforming  ^ds.  However,  the  strategy  (2)  is 
preferable  since,  in  this  case,  the  viscous  solution  is  independent  of  the  inviscid  one. 
Calculations  were  made  using  strategy  (2)  and  a  fixed  grid,  but  were  found  to  be 
computationally  unstable.  Also,  it  was  found  that  the  wave  field  violated  the  radiation 
condition.  Presently,  an  unsteady  approach  is  under  development  in  which  the  solution  is 
initiated  from  a  condition  of  zero  freestream  velocity  and  accelerated  to  a  steady-state 
nonzero  value.  A  laminar-flow  code  is  under  development  for  the  initial  calculations. 

(c)  Extensions  are  being  made  of  the  interactive  approach  for  yawed-ship 
calculations.  Figures  5  and  6  provide  SPLASH  results  for  the  Series  60  Qb  =  .6  ship 
model  for  Fr  =  0  and  .316,  respectively,  and  yaw  angles  y  =  0,  5,  and  10  deg.  Both 
barebody  and  Fr  =  0  displacement  body  results  are  shown.  Presently,  the  results  are  being 
evaluated  through  comparisons  with  the  experimental  data. 

Extensions  are  also  in  progress  of  the  viscous-flow  code  for  interactive  nonzero  Fr 
yawed-ship  calculations.  Initially,  zero  Fr  calculations  are  being  performed  to  validate  the 
extensions.  Large-domain  calculations  have  been  completed  and  are  considered 
satisfactory  (Van  and  Stem,  1991).  Presently,  extensions  are  being  made  to  resolve  the 
bow  flow  in  conjunction  with  interactive  calculations,  first  for  Fr  =  0,  and  then  for  nonzero 
Fr. 


4.  Interactions  with  Navy  Agencies 

1.  In  response  to  I>r.  Rood's  letter  of  20  November  1991,  an  abstract  has  been 
submitted  for  a  proposed  presentation  at  the  1992  Free-Surface  Vorticity  Workshop 
entitled  "Solid-Huid  Juncture  Boundary  Layer  and  Wake  with  Waves." 

5.  Publications 

Longo,  J.,  F.  Stem,  and  Y.  Toda,  (1991),  "Mean-Flow  Measurements  in  the  Boundary 
Layer  and  Wake  and  Wave  Field  of  a  Series  60  Cb  =  -6  Ship  Model  -  Part  2:  Scale-Effects 
on  Near-Field  Wave  Patterns  and  Comparisons  with  Inviscid  Theory,"  submitted  to  J.  of 

Ship  Rcsgargh- 

Stem,  F.,  Choi,  J.  E.,  and  Hwang,  W.S.,  (1991),  "Effects  of  Waves  on  the  Wake  of  a 
Surface-Piercing  Flat  Plate:  Experiment  and  Theory,"  Iowa  Institute  of  Hydraulic 
Research,  The  University  of  Iowa,  IIHR  Report  No.  353;  also  submitted  to  J.  of  Ship 
Research. 

Tahara,  Y.  and  Stem,  F.,  (1992)," Validation  of  an  Interactive  Approach  for  Calculating 
Ship  Boundary  Layers  and  Wakes  for  Nonzero  Froude  Number,"  in  preparation  for 
submission  to  Computers  and  Fluids. 

Toda,  Y.,  Stem,  F.,  and  Longo,  J.,  (1991),  "Mean-Flow  Measurements  in  the  Boundary 
Layer  and  Wake  and  Wave  Field  of  a  Series  60  Cb  =  -6  Ship  Model  for  Froude  Numbers 
.16  and  .316,"  Iowa  Institute  of  Hydraulic  Research,  The  University  of  Iowa,  lEHR  Report 
No.  352;  also  to  appear  in  Proc.  of  The  Second  Osaka  International  Colloquium  on 
Viscous  Fluid  Dynamics:  also  submitted  to  J.  Ship  Research  under  title  "Mean-Flow 
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Measurements  in  the  Boundary  Layer  and  Wake  and  Wave  Field  of  a  Series  60  Cb  =  .6 
Ship  Model  -  Part  1:  Froude  Numbers  .16  and  .316." 

Van,  S.-H.  and  Stem,  F.,  (1991),  "Computation  of  Viscous  Flows  around  Wigley  Hull 
Advancing  with  an  Angle,  of  Attack,"  presented  at  the  fall  meeting  SNAK. 

L _ Miscellaneous 

1.  Dr.  Stem  gave  an  invited  presentation  to  the  Paitnership  for  Americas  Cup 
Technology  on  16  September  entitled  "Wav^oundary-Layer  and  Wake 
Interaction." 

2.  Dr.  Stem  attended  and  gave  a  presentation  (Toda  et  al.,  1991)  at  The  Second 
Osaka  International  Colloquium  on  Viscous  Ruid  Dynamics  in  Osaka,  Japan  on  27 
through  30  September  1991.  During  this  period  Dr.  Stem  held  meetings  with  Dr. 
Toda  concerning  the  Osaka/IIHR  cooperative  project  on  propeller-hull  interaction. 

3.  Dr.  Stem  visited,  by  invitation,  CSSRC  (Wuxi,  China)  and  MARIC  (Shanghai, 
China)  and  gave  presentations  entitled  "Computational  Ruid  Dynamics," 
"Wave/Boundary-Layer  and  Wake  Interaction,"  and  "Viscous  Propulsor 
Hydrodynamics." 

4.  A  proposal  was  submitted  for  300  hours  of  supercomputer  time  on  the  NAS 
system  for  FY  92  in  support  of  future  related  research. 


Respectfully  submitted  by. 


cc:  Mr.  Robert  L.  Silverman,  ONR,  Seattle 


Fred  Stem 

Principal  Investigator 
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1  zero  gradient  at  z  =  0 

2  pressure  w/o  curvature  effect  at  z=0 

3  stress  w/o  curvature  effect  at  z=0 

4  exact  at  z  =  0 

5  pressure  w/o  curvature  effect  at  z=r\ 

6  stress  w/o  curvature  effect  at  z=T| 

7  exact  at  z  =  Tj 
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Free-surface  conforming  grid : 

Mean  free-surface  grid : 


Figure  2  Comparison  of  the  Relative  Error  for  Various  treatments 
of  the  Free-Surface  Boundary  Conditions 


(b)  vorticity  difference 

Figure  3  Difference  Between  Zero  and  Nonzero  Steepness  Solutions 


(c)  depthwise  vorticity  equation 


Figure  4  Vorticity  Equation 


PRESSURE  CONTOURS  A=  O.OUO 


Figure  5  SPLASH  Results  for  0,  5,  and  10  Degree  Yaw  :  Fr 
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•  FREESl^RFACE  CONTOI  RS  ^=  0.0000  FREESURFACE  CONTOURS  i=  O.OOOT  FREESURFACE  CONTOURS  ^  0.0015 


Figure  6  SPLASH  Results  for  0, 5,  and  10  Degree  Yaw :  Fr  -  .316 


